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Spitzer photometry and spectroscopy of the star HD69830 reveal an excess 
of emission relative to the stellar photosphere between 8 and 35 /im dominated 
by strong features attributable to crystalline silicates with an emitting surface 
area more than 1,000 times that of our zodiacal cloud. The spectrum closely 
resembles that of the comet C/1995 01 (Hale-Bopp). Since no excess is detected 
at 70 /im, the emitting material must be quite warm, be confined within a few 
AU of the star, and originate in grains with low long-wavelength emissivity, i.e. 
grains much smaller than 70 /zm/27r ~ 10 /im. The strong mineralogical features 
are evidence for even smaller, possibly sub-micron sized grains. This small grain 
size is in direct contrast to the 10-100 /xm grains that dominate the relatively 
featureless spectra of our zodiacal dust cloud and most other main sequence stars 
with excesses. The upper limit at 70 /im also implies that any Kuiper Belt analog 
must either be very cold or be less massive than ~ 5 times our own Kuiper Belt. 

With collisional and Poynting-Robertson drag times of about a thousand 
years for small grains, the emitting material must either: a) be created through 
continual grinding down of material in a dense asteroid belt; or b) originate in 
cometary debris arising from either a single "super-comet" or a very large number 
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of individual comets arriving from a distant reservoir. In the case of a cometary 
origin for the emission, the mass requirements for continuous generation by many 
individual comets are unreasonable and we favor the capture of a single "super 
comet" into a 0.5-1 AU orbit where it can evolve a large number of small grains 
over a 2 Myr period. 

Subject headings: (stars:) circumstellar matter, debris disks, Kuiper Belt, comets, 
infrared, extra-solar planets 

1. Introduction 

The debris disk phenomenon, first discovered by the Infrared Astronomy Satellite (IRAS) 
(Aumann et al. 1984; Gillett 1986; Backman & Paresce 1993) is of great interest because of 
the clues these disks offer toward understanding the formation and evolution of planetary 
systems. Toward this end, aspects of debris disks have been studied with the Infrared Space 
Observatory (ISO) (Habing et al. 2001; Decin et al. 2003), optical imaging (e.g., Smith & 
Terrile 1984; Heap et al. 2000), and in the submillimeter (Holland et al. 1998; Greaves et 
al. 1998; Dent et al. 2000). 

Information on the nature of the dust grains in these disks is particularly impor- 
tant. Spectroscopic observations toward luminous Herbig Ae/Be stars such as HD100546 
(Grady et al. 1997; Malfait et al. 1998; Meeus et al. 2001) and in young debris stars such 
as (3 Pictoris (Telesco and Knacke 1991) and 51 Oph (Fajardo-Acosta et al. 1993) revealed 
dust composed, at least in part, of small (sub-micron) grains of crystalline silicates such 
as forsterite and enstatite. The similarity of these spectral features to those seen in comet 
C/1995 01 (Hale-Bopp) (Crovisier et al. 1996; Griin et al. 2001; Wooden et al. 2000) sug- 
gests that this circumstellar material may represent cometary debris. Yet these spectral 
features are not present in all or even most debris disks. More than a dozen classic debris 
disks (including Fomalhaut) examined by Spitzer (Jura et al. 2004; Stapelfeldt et al. 2004) 
show little or no spectral structure, suggesting that the grains in these systems are larger than 
~10 fim. These grains may be similar to those in our own zodiacal cloud which are predomi- 
nantly larger than 10-100 /zm with only a small admixture of smaller silicate grains yielding a 
weak 10 /xm emission feature (Kelsall et al. 1998; Fixsen & Dwek 2002; Reach et al. 2003). 

In this paper we report the detection of a large excess due to hot grains orbiting the star 
HD69830. These grains are dominated by a population of crystalline silicates with prominent 
features in the IRS spectrum. 
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2. Properties of HD69830 

HD69830 (Gliese 302, BD-12 2449, HR 3259) is a bright (V=5.95 mag), nearby (12.6 
pc), main sequence dwarf of solar metallicity, [Fe/H] =0.015 (Cayrel de Strobel et al. 2001). 
Its spectral type has been variously classified as between G8 (Skiff 2003) and K0 (Song et 
al. (2000) and other references in SIMBAD). We adopt a K0V spectral type throughout 
this paper which implies an effective temperature of 5150 K, a mass of 0.8 M , and a total 
luminosity of 0.45 L (Cox 2000). Using isochrone fitting, Li absorption, Ca II H&K line 
activity, weak X-ray emission, and space motions, Song et al. (2000) suggested an age older 
than 0.6 Gyr and up to 2 Gyr with a preference for the older age. Examination of Ca II 
H&K activity indices leads Marcy and Fischer (private communication, 2005) and Wright 
et al. (2004) to prefer an age of 2-4.7 Gyr while IUE spectra suggest an age of 3 Gyr 
(Hufnagel & Smith 1994). We adopt 2 Gyr as the probable age for this star. 

No planets have been detected orbiting HD69830 using the radial velocity technique 
with a limit on any planet with a mass, M sin(i) > 0.5 Mj up located within a distance of 
3 AU of the star (Marcy and Fischer, private communication). Although Mannings and 
Barlow (1998) noted the presence of a small "hot" excess based on IRAS observations at 
25 /im (3.8<r, Table 1), IRAS did not detect the star at either 60 or 100 /an. The star was 
not observed with ISO. 



3. Spitzer Observations 

3.1. Observations and Data Reduction 

Spectra of HD69830 were obtained with Spitzer on 2004 April 18 (AOR 4016640). These 
measurements were part of a program that has observed 36 main-sequence stars to look for 
disks (Beichman et al. 2005; Bryden et al. 2005) using the Infrared Spectrograph (IRS) 
(Houck et al. 2004). Standard IRS Staring mode observations were made with the Short 
Low Order 1 (SL1; 7-14 /im), Long Low Order 2 (LL2; 14-20 pm), and Long Low Order 1 
(LL1; 20-35 /im) modules. Each star was observed at two positions along the slit, called Nodi 
and Nod2. Data were processed by the Spitzer Science Center (SSC) to produce calibrated 
images of the spectrometer focal plane. The SSC SPICE software was used to extract spectra 
from these images. Because of the need for careful subtraction of the stellar continuum to 
detect a faint excess, additional steps were taken in producing the final spectra: 

• To ensure proper sky subtraction, differences were formed between the two Nod posi- 
tions to produce two sky-subtracted images (Diff 1 =Nod 1 — Nod 2 and Diff 2 =Nod 2 — Nodi). 
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SPICE was used to extract 1-D spectra from each difference image using default set- 
tings and calibration. 

• The first few points at the beginning and end of the spectrum from each module were 
typically unreliable, as were a few bad pixels flagged in the SSC processing. These 
were rejected. This effect was particularly noticeable at the long wavelength end of 
SL1 (A > 14 /mi). Since the short wavelength end of LL2 overlaps SL1, there is no gap 
in the final spectrum. 

• Longward of 14 /zm, photospheres of solar type stars are smooth and do not differ 
greatly from a Rayleigh- Jeans blackbody (Kurucz 1992; Castelli & Kurucz 2003). We 
used this fact to improve the pixel-to-pixel calibration in the extracted SSC/SPICE 
spectra. 35 of the 36 stars observed with IRS as part of this program showed no 
evidence for an excess at either 24 or 70 /im with MIPS, but did show consistent 
pixel-to-pixel deviations from a smooth photospheric model. From these stars we 
formed ratios of the extracted IRS spectra to the Kurucz model appropriate for the 
effective temperature and metallicity of each star. From the average of these ratios at 
each wavelength, we created a "super-flat" response curve for the SL1, LL2 and LL1 
modules. Because of possible variability in the flat response in the SL1 module, we 
used only the 8 other stars observed during the same instrument campaign (April 2004) 
as HD69830 to generate the super-flat for this module. For the LL modules we used 
data from all 35 stars. HD69830 was excluded from the creation of the super-flats. 

The corrections derived from the comparison of extracted spectra and the photospheric 
models for the 3 modules are small: 

• Shortward of 14 /im, the SL1 super-flat typically has values in the range 0.97-1.03, or, 
equivalently, deviations from unity response of ~3%; a few pixels have values deviating 
from unity by 10%. The values at each pixel of the SL1 super-flat have a dispersion 
around the pixel response (averaged over all 8 stars) of (J pop r ^ J 2% and (7 mean ^ 

0.4%. In 

the LL1 and LL2 modules, the typical super-flat values are in the range 0.98-1.02 with 
dispersions (averaged over all 35 stars) of a pop ~ 1 — 2% and a mean ~ 0.2 — 0.3%. We 
divided the extracted spectra by the appropriate "super-flat" to remove any residual 
calibration variations. The dispersion in the super-flat {cr pop ) is an indication of the 
limiting systematic noise in removing the photospheric contribution to the signal from 
these stars. 

• It was necessary to adjust the absolute flux scale of all of the IRS/SSC spectra by a 
factor of 1.20±0.01 to fit the photospheric models obtained by fits to shorter wavelength 
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(0.4-2.2 /mi) data. With this adjustment the SL1 data are consistent with the IRAS 
12 /iin data and the LL2 data are consistent with the IRAS/MIPS data at 24 /im. The 
origin of this absolute calibration factor, e.g. pointing errors or an imperfect aperture 
correction, is being pursued but is not relevant to the present investigation. 

3.2. Results 

The overall Spectral Energy Distribution (SED) of HD 69830 including photometry 
(Table 1) and the IRS spectra is shown in Figure 1. A prominent excess is obvious in all 
three IRS modules reaching 60% over the photospheric continuum at 35 /mi. A number 
of broad spectral features, discussed below, are prominent in the data. Figure 2 compares 
the fractional deviation from a smooth photosphere for HD69830 with the spectra of two 
stars, HD 173667 (AOR 4024320) and HD142373 (AOR 4023040), that show no excesses at 
either IRAS or MIPS wavelengths. As is clear from this figure, the size and spectral content 
of HD60930's excess are large compared to the uncertainties and will not be qualitatively 
altered by further refinements in the spectral extraction. Figure 2 can also be taken as a 
limit to the deviations which might be found toward any star, approximately 2-3%. 

Included in Figure 1 and presented in Table 1 is MIPS photometry (AOR 4041728) 
as well as color-corrected IRAS data from the Faint Source Catalog (Moshir et al. 1990). 
The MIPS data were reduced using the MIPS/GTO DAT software (Gordon et al. 2005) and 
show a prominent (50%) excess at 24 /im, confirming the excess hinted at in the IRAS data 
at 25 /im (Mannings & Barlow 1998). There is no evidence for extended emission in the 
24 /im image. The consistency of the IRAS and Spitzer/MIPS data at 24 and 25 /mi suggests 
that there has been no strong variation on the 20 year timescale. 

As discussed elsewhere, HD69830 is the only star out of a sample of 84 showing a 
significant 24 /mi excess (Beichman et al. 2005; Bryden et al. 2005). The uniqueness of 
HD69830 in this sample, 1 of 84 or 1.2%, confirms the rarity of hot excesses around ma- 
ture FGK stars first observed by IRAS (Fajardo-Acosta, Beichman & Cutri 2000) and ISO 
(Laureijs et al. 2002). But HD69830 is even more exceptional in that its 24 /im excess 
is not accompanied by a longer wavelength excess as is normally the case for debris disk 
sources (Mannings & Barlow 1998; Backman 1998). The measured value of the HD69830 
excess at 70 /im is 1 ± 3 mJy compared to the predicted photospheric level of 18 mJy, or 
F ^cess I F star <50% (3 times the 1 a statistical uncertainty combined in quadrature with a 
20% overall calibration uncertainty) at 70 /mi. The lack of a 70 /im excess is evidence for 
the presence of small grains with low, far-infrared emissivity. 
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The upper panel of Figure 3 shows the IRS data with the photospheric component 
subtracted to reveal the details of excess. Uncertainties at each wavelength were derived 
using the largest of two terms: a) the quoted errors in the individual SPICE extractions; and 
b) the difference between the corrected spectra obtained at the two different nod positions. 
To give the total uncertainty, the preceding value is combined in quadrature with a systematic 
term equal to 1% of the photospheric emission to account for errors in the subtraction of the 
photosphere. The first two terms reflect statistical uncertainties in the spectrum and become 
worse at the longest wavelengths where the signal to noise is poorest. The third term, being 
proportional to the stellar brightness, decreases in importance at longer wavelengths. 

4. The Excess Toward HD69830 

The lower panel of Figure 3 shows a spectrum obtained by the Infrared Space Observa- 
tory (ISO) of the comet Hale-Bopp (Crovisier et al. 1996) and scaled as described in §4.2 to 
match the overall slope of the Spitzer data. The detailed similarity between the two spectra 
is remarkable and will be used as a guide to the identification of features in the spectrum of 
HD69830. 



4.1. Features in the Spectrum of the excess 

Table 2 identifies some of the features in the excess. Most prominent are broad fea- 
tures at 9.4-12 //in, 16 /im, 19.3 /im and 23.8 /im. Other features of moderate strength 
include a plateau at 10 /an, a plateau at 27 /im, and a weak feature of low signal-to-noise 
ratio around 33-34 /im. As noted in Table 2 and Figure 3, small Mg-rich crystalline silicate 
grains 1 , e.g. enstatite and forsterite, have spectral features that match those in the spectrum 
remarkably well (Jager et al. 1998; Koike, Shibai, & Tuchiyama 1993). These are seen in 
the ISO and ground based spectra of Hale-Bopp (Crovisier et al. 1996; Crovisier et al. 1997; 
Wooden et al. 1999; Wooden et al. 2000). While not dominant in the spectrum of the zodi- 
acal cloud, these crystalline silicates are found in interplanetary dust particles as well as in 
meteorite inclusions (Yoneda et al. 1993; Bradley 2003). 

A detailed examination of the spectrum reveals that forsterite matches the 9-11 /im 



1 As discussed in Wooden et al. (1999 and 2000) olivines are minerals with [Mg y Fei_ y ]2SiC>4 and < y < 1; 
forsterite is an Mg-pure olivine with y=0. Pyroxenes are minerals with [Mg^Fei^JSiOa and < x < 1; 
enstatite is a Mg-pure pyroxene with x=0. 
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feature, while enstatite more closely matches the ratio of the 19 jum and 23.8 /im features 
(Henning 1998). Features due to crystalline olivines are far more prominent than those due to 
pyroxenes, although an inflection point around 27 /im might be due to crystalline pyroxene. 
The 9.3 /im feature characteristic of small, amorphous silicate grains that is prominent in 
Hale-Bopp and seen weakly in our own zodiacal cloud (Reach et al. 2003) is not seen at all 
in HD69830. Since the emission between the features does not drop to zero, there must be 
some source of continuum opacity, possibly larger grains, in both HD69830 and Hale-Bopp. 
In summary, we attribute the excess toward HD69830 to a mixture of grains similar to those 
seen in Hale-Bopp and very different from those responsible for the relatively featureless 
spectra of typical debris disks. 



4.2. Model for the Excess Emission 

The simplest model for the emission (Figure 3) scales the ISO spectrum of Hale-Bopp 
by the ratio of two blackbodies: 

F U (H D69830) oc F u (HaleBopp) x B u (A00K) / B u (T(Hale Bopp) ) (1) 

where B U (T) is the Planck function and T (Hale-Bopp) =207 K is the equilibrium temperature 
of a small crystalline silicate grain 2.9 AU from the Sun. The value of 400 K is fitted by eye 
to provide a reasonable fit to the HD69830 excess and corresponds to an equilibrium distance 
of ~0.5 AU for small grains. Although the single temperature scaling gives a remarkably 
good representation of the data, we next develop a more realistic model that constrains the 
location, size, and composition of the emitting material. 

The presence of pronounced spectral features on top of a weak continuum suggests 
optically thin emission dominated by dust grains comparable to or smaller than the wave- 
length: a < X c /2tt ~ 1 /im where A c ~ 8 /im is the shortest wavelength showing spectral 
features. The close analogy to Hale-Bopp requires the existence of grains as small as 0.25 /im 
(Wooden et al. 1999; Wooden et al. 2000). While fluffy, low density grains with a size around 
1 /im have been suggested to fit the spectrum of Hale-Bopp (Li & Greenberg 1998), for this 
preliminary analysis we adopt the following empirical model: 

• Synthesize empirical absorption cross sections, Q a b s , derived from the scaled Hale- 
Bopp observations (Wooden et al. 1999; Wooden et al. 2000), Q a bs oc F v /B v (T H b), 
with Tub = 207 K and scaled to match laboratory measurements of crystalline olivine 
at 11.3 /im (Mukai & Koike 1990). 
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• Subtract optical constants for 0.25 /iin amorphous silicate grains (Draine & Lee 1984; 
Weingartner & Draine 2001) from the optical constants derived from Hale-Bopp since 
the 9.3 /im feature produced by this type of grain is not seen in HD69830. Examination 
of Hale-Bopp emission models, e.g. Figure 6 in Wooden et al. (1999) suggests that 
amorphous grains are at least a factor of five less abundant in HD69830 than in Hale- 
Bopp. 

• Synthesize empirical absorption cross sections at 70 fim as described above using the 70 
fim data for Hale-Bopp from Griin et al. (2001). The assumption that all the emission 
comes from grains following the same equilibrium temperature law is naive but gives 
a self-consistent estimate of the 70 /im emission. 

• Assuming a flat disk geometry for convenience, calculate the temperature gradient 
through radiative equilibrium of 0.25 fim grains orbiting a K0 star and integrate the 
surface brightness profile, I v {r) = r(r)B v (T(r)) over a region extending from the grain 
sublimation distance (at which T~1500 K and corresponding to r ~ 0.06 AU) out to 
an outer radius r max . 

• Let the radial dependence of the optical depth follow a power law distribution, r(r) oc 
r a . For comparison, a = — 0.4 in our solar system. 

A x 2 search was made through a broad range of values for the free parameters in the 
model: the surface area of the small grains derived from the Hale-Bopp spectrum, the power 
law index of the surface density profile, and r max . The parameters of the best fit model 
are given in Table 3. As shown in Figure 4, the model reproduces many, but not all, of the 
details in the IRS spectrum. Some of problems include a shoulder at 10 fim that cannot be fit 
without throwing off the fit at 11 fim; an excess between 13-15 fim not seen in Hale-Bopp (the 
feature does not correspond to known Mg- or Mg/Fe-rich silicates, but resembles unidentified 
features in the spectrum of R Cas (? )); subtle variations around 19 fim and 24 fim where 
4 or 5 separate bands of forsterite and enstatite make up a broad feature sensitive to small 
composition and/or temperature differences (Bowey et al. 2001). Despite these small but 
significant differences in detail, the Hale-Bopp based model provides an excellent analogue 
for the HD69830 debris disk. 

In the final fit (Table 3 and Figure 4) we constrained the power law index of the radial 
surface density power law to have a value of —0.4 since varying this parameter away from 
-0.4 did not significantly improve the overall results. Models with steeper power law indices 
and r max > 5 AU were consistent with the limit to the 70 fim excess, but had poor fits 
in the 7 — 35 fim region with x 2 > 10- Models with a = —0.4 and r rnax > 1 AU also 
had poor x 2 values and were inconsistent with the 70 fim limit. Although grain size was 
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not an explicit free parameter in the model because we adopted opacities directly from the 
Hale-Bopp spectrum, we found that grains larger than about 5 pm were inconsistent with 
the limits on the 70 pm excess for any reasonable combination of model parameters. 

The fractional luminosity in the excess, integrated between 7 and 35 pm is Ldust/L* = 
2 x 1CT 4 , making this system a factor of two brighter than Fomalhaut (a PsA; Backman and 
Paresce 1993; Stapelfeldt et al. 2004) in terms of fractional dust luminosity and more than 
1,000 times brighter than our own zodiacal cloud. The surface area of emitting material in 
the HD69830 cloud is 1.3 x 10 3 times that of our own zodiacal cloud, ~ 2 x 10 20 cm 2 within 
2-3 AU (Backman 1998). In terms of a fractional surface area of emitting material, our solar 
zodiacal cloud has £/ rac ~ 1 x 10—7 which compares to the fractional surface area of 0.25 
pm grains within 1 AU of S /mc ~ 3 x 10" 4 for HD69830. 

If the emitting grains all have radius a = 0.25 pm with density p = 3 g cm~ 3 then the 
mass of the emitting material (Table 3) would be | Spa = 2.7 x 10 19 g or 4.6 x 10~ 9 M® 
where S = N tot ira 2 is the emitting surface area of N tot grains. However, since HD69830 is 
probably surrounded by a broad distribution of grain sizes, the total mass could be much 
larger. For a grain size distribution n(a) = n (a/a )~ 3 ' 5 , corresponding to predictions for 
a collisional cascade (Dohnanyi 1969), the ratio of total mass, M, to surface area, S, is 
(M/E) ~ ^p^a-minO-max where a m i n is the minimum grain size and a max is the maximum 
grain size. For a m i n = 0.25 pm and a max = 10 km, the total mass of material associated with 
the disk emission (Table 3) is 5.4 x 10 24 g or 0.9 x 10~ 3 M e ; given the strong dependence 
on the unknown grain size distribution, this extrapolation is obviously very uncertain. This 
number should be compared to the cumulative mass of asteroids up to 10 km in our solar 
system, 0.014 x 10~ 3 M ffi (Bidstrup et al. 2004; Krasinsky et al. 2002), indicating a factor of 
64 times more mass in a putative "asteroid" belt around HD 69830 than in our own asteroid 
belt (up to 10 km size). We analyze this estimate in the context of grain production and 
loss mechanisms in §4.4. 

4.3. Limits on a Kuiper Belt 

One of the most striking aspects of HD69830 is the lack of a long-wavelength ex- 
cess to accompany the 8-35 pm excess. Treating the upper limit at 70 pm just like the 
upper limits toward other main sequence stars observed by MIPS (Beichman et al. 2005; 
Bryden et al. 2005) yields an estimate for the amount of cold dust. The peak of the cold 
dust emission comes at 70 pm for material at a temperature of 50 K located at ~ 50 AU 
from the star (for small grains; perhaps a factor of 2 closer in for larger grains). In this 
case, L dust /L* < 10~ 5 (5600 K/T*) 3 (F dust /F*) ~ 5 x 10~ 6 (3d). This is only 5 times larger 
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than the upper end of the range of values, 10 -6 to 10~ 7 , predicted for the Kuiper Belt region 
beyond 30 AU in our solar system (Stern 1996) and less than the IRAS constraint, < 10~ 5 
(Aumann & Good 1990). The HD69830 limit on cold material is uncertain as it is based 
on only a single wavelength; observations in the submillimeter would further constrain the 
potential reservoir of distant cold material (Greaves et al. 2004). 



4.4. Origin of the Emitting Grains 

Nothing in the properties of HD69830 suggests that this star is anything other than a ma- 
ture, main sequence dwarf. Thus, there is no reason to attribute its excess to dust associated 
with the weak mass loss often seen toward M and K giants (Beichman et al. 1990). Spectral 
features in the excess emission around such stars have been identified with small crystalline 
grains condensing in the quiescent, outflowing gas (Henning 1998; Dijkstra et al. 2003), but 
since HD69830 is a quintessential dwarf star, we do not consider this possibility further. 
In the absence of any purely stellar phenomenon, we attribute the excess to orbiting dust 
derived from the solid material left over from the formation of the star itself, i.e. a debris 
disk. We now investigate the source of this dust in more detail. 



4-4- 1- Timescales 

The rarity of main sequence stars with hot excesses and without an associated colder 
debris cloud (<1%) suggests that the conditions for this excess are found around fewer than 
1 star in 100 or that it is a transitory (or low duty cycle) event lasting no more than 1% of 
the 2 Gyr lifetime of HD69830. 

What are the natural timescales of the debris clouds in HD69830, and, for comparison, in 
our solar system? A debris cloud continually loses material due to collisional grinding down 
of grains to a size small enough that radiation pressure and/or Poynting Robertson (PR) 
drag ultimately remove them from the system. The relevant timescales (Burns et al. 1979; 
Backman & Paresce 1993) are 



torbit 

82/rac 8Ej rac il£, 



Tcoll = ^ = 77v-i 77oI y r ( 2 ) 



where £/ mc is the fractional surface area of the dust (Backman & Paresce 1993) and the PR 
drag time (Backman & Paresce 1993) 

tpr = 700 ag rain p r 2 AU ^ yr (3) 

Li* 
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where p (g cm 3 ) is the particle density and the grain size is in /zm. 

Table 4 compares relevant parameters for our solar system's zodiacal cloud and that 
of HD69830. In our solar system the critical (shortest) timescale is the PR time, r PR = 
2 - 6 x 10 5 yr for 10-30 /im grains (Fixsen & Dwek 2002; Reach et al. 2003). In HD69830 
the collisional timescale is 400 yr and the PR timescale about 700 yr for 0.25 /im grains at 
1 AU. The grains will grind themselves down through collisions until they are small enough 
to be lost through PR drag, i.e. until Tp R ~ t co u. Following Sheret et al. (2004; see their 
Figure 3) we note that radiation blowout does not operate efficiently for low luminosity 
stars like HD69830 and e Eri, particularly for non-idealized grains that are either porous 
or non-spherical (Li & Greenberg 1998). Collisions followed by loss via PR drag are likely 
to dominate the grain distribution around HD69830. Although the uncertainties on the 
HD69830 timescales are large (at least a factor of two due to unknown grain size and orbital 
location), it is clear that small grains must be replenished on roughly the < 1,000 yr timescale 
or the emission associated them will disappear. 



4 -4 -2. An Enhanced Zodiacal Dust Cloud 

The probable source for the small grains seen toward HD69830 is a debris cloud located 
within 1 AU of the star that is replenished on timescales of 10 3 yr. For comparison, our zodi- 
acal cloud consists of relatively large grains (10-100 fxm) extending to 3 AU and maintained 
in rough equilibrium, with occasional upward spikes, through asteroid collisions every few 
million years (Fixsen & Dwek 2002; Dermott et al. 2002; Dermott et al. 2002). 

The mass in 0.25 /im dust grains orbiting HD 69830 (§4.2), 4.6 x 10~ 9 M e , is destroyed 
in T co u ~ t pr ~ 0.5 — 1 x 10 3 yr for a mass loss rate of 0.5 — 1 x 10 -11 M e yr -1 . This loss 
mechanism represents the end point for the small grain debris generated by the collisional 
cascade of larger bodies into smaller ones. If this mechanism operates continuously over 2 
Gyr, then the total mass of solid material lost to the system over 2 Gyr would be > 10~ 2 
M e , or more than 500 times the amount of material in own asteroid belt. Note that the 
estimate of solid material lost over time is also considerably greater than the present day 
mass estimate, 1 to3 x 10~ 3 M e inferred from integrating an a -3 ' 5 distribution up to either 
10 or 100 km bodies (§4.2). Either we are seeing the decaying remnant of a more massive 
disk or the mass loss is episodic and the average mass loss rate is not as great as present day 
observations would suggest. 

Another way to estimate the amount of material present is to examine the dust pro- 
duction rate. If, as noted above (Sheret et al. 2004), radiation blowout is not effective 
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around HD69830, then PR drag is the dominant loss for dust in HD69830 and the number of 
emitting grains in the disk, n gr , is simply proportional to the square of the number of large 
bodies (comets or asteroids), iV c , that initiate the collisional cascade leading to these small 
grains: n gr oc N% (Dominik & Decin 2004). Presently, our zodiacal cloud has a fractional 
surface density around 10~ 7 , but immediately after a major collision, the cloud would have 
a larger population of small grains given by the n(a)oc a~ 3 - 5 law appropriate to a collisional 
cascade (Dohnanyi 1969). The surface area of emitting material and the resultant optical 
depth of such a cloud is oc J a 2 a -3 - 5 da oc a -0 5 . Extending the small grain cutoff from its 
present 10 fim to 0.25 /im would temporally increase the emitting surface area by a factor of 
at least a factor ~ 6 and probably more since a collision would increase the amount of dust 
of all sizes. If we compare the cloud around HD69830 with our solar system (enhanced to in- 
clude short-lived 0.25 /iin grains for this calculation), we need Sj mc (HD69830) / Sj mc (Solar 
System) = 3 x 10~ 4 /6 x 10~ 7 = 500 times more small grains orbiting HD69830 than our 
zodiacal cloud would have immediately after a major collision. 

Since, as argued above, S/ rac oc N% (Dominik & Decin 2004), the number of large 
bodies required to replenish the dust would be N c oc V500 = 22 more than in our solar 
system. By this reasoning the total mass of a putative asteroid belt would be 0.3 xl0~ 3 
M e which compares favorably with the value of 1 xl0~ 3 M e derived from integrating an 
a -3.5 g j ze distribution up to 10 km. Although both mass estimates would increase if the 
upper size mass limit were increased to 100 km, the values are considerably smaller than 
the 10~ 2 M e estimated by integrating the instantaneous mass loss rate over the maximum 
10 7 yr lifetime of the excess. The simplest hypothesis is that the numbers are consistent 
within the uncertainties (particularly since the derived masses depend on the unconstrained 
upper end of the integral over the particle size distribution, 10 ~ 100 km). If the difference 
between the instantaneous and long term mass rates is significant, then we could be seeing 
a transient peak in surface area, possibly from a very recent collision (< tpr ~ 10 3 ) yr. 

The exact details of an asteroid belt model will have to be established by a combination 
of improved modeling to constrain grain properties and location more precisely, as well as 
more realistic modeling of the dynamics. In any event, HD69830 could have an asteroid 
system with much higher density located much closer to the star than in our solar system. 
Pursuing the analogy with our solar system suggests that an as-yet-unseen planet located 
outside r max could truncate the disk and confine most of the dust to radii interior to the 
planet's orbit (Menou & Tabachnik 2003). Such a planet could also enhance the lifetime of 
dust by trapping it into mean motion resonances (Kuchner & Holman 2003). 
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4-4-3- A Swarm of Comets as a Source of Small Grains 



An alternative to emission from small grains in an enhanced asteroid belt is emission 
from cometary ejecta. Levison et al. (2001) discuss possible mechanisms for stirring up the 
Kuiper Belt including scattering of material into the inner solar system by the late formation 
of distant planets or the destabilization of Kuiper Belt objects due to sudden changes in the 
orbital parameters of one or more giant planets. Alternatively, some perturbing event might 
send a single large object located far from the star, i.e. a 1,000 km sized object like Sedna 
(Brown et al. 2003), into the inner HD69830 system where it becomes a "super-comet". We 
examine these two possibilities for HD69830 below. 

What would it take to reproduce the emission from HD69830 using cometary debris? 
Consider the F^(11.3;um) = 7300 Jy flux density of Hale-Bopp measured in a 9.3" beam 
when Hale-Bopp was located at = 1.1 AU from the Sun and Ahb = 1-6 AU from the 
Earth (Williams et al. 1997). If we scale this to Hale-Bopp's perihelion values of = 0.9 
AU and A HB = 1.5 AU according to r^ 4 A~ 1 (Gehrz & Ney 1992; Wooden et al. 1999) and 
to a radius of 20" (Hayward, Hanner, & Sekanina 2000), then the total 11.3 emission 
of this object would be ~ 6 x 10 4 Jy. Placed at the A(HD69830) = 12.6 pc distance of 
HD69830, it would take 



comets identical to Hale-Bopp as it appeared at 0.9 AU to reproduce the 0.013 Jy flux density 
measured by Spitzer at 11.3 /iin. Reducing the periastron distance to 0.5 AU consistent with 
the ~ 400 K temperature determined in the modeling would increase the brightness of the 
cometary material and thus decrease their required number to ~ 5 x 10 5 Hale-Bopps. 

Consider a steady stream of comets on Period = 10 3 yr orbits arriving at the inner 
HD69830 system at a rate of rj comets yr -1 . If each comet is bright while it is within 3 AU 
for a time t 3AU ~ 1 yr and makes N Pass ~ 100 active passes through the inner HD69830 
system (Griin et al. 2001; Wooden et al. 2000), then in steady state N H b = V^au comets 
contribute to the observed excess at any one time; N ac u ve = r]Period contribute to the excess 
whenever they come close to the star; and N tota i = 77 Lifetime /Np ass comets contribute to 
the excess during the putative 10 7 year lifetime of the excess. If, as argued above, we need 
5 x 10 5 Hale-Bopps at any instant, then i] = 5 x 10 5 comets yr" 1 . With these numbers, 
Nactive = 5 x 10 8 and N to tai = 5 x 10 10 over the duration of this excess. With a total mass for 
Hale-Bopp of ~ 1 x 10 20 g, assuming a density of 1.5 g cm" 3 for a mixture of grain and ice 
and a radius of ~ 25 km (Weaver & Lamy 1997), then the total mass in comets entering the 
inner HD69830 system during this 10 7 yr period would be ~ 900M e , unreasonably large for 
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any residual Kuiper Belt. Even if more normal, smaller comets were adopted, e.g. a factor 
of 10 smaller than Hale Bopp, the total mass would still require ~ 90M e in comets which 
still seems large given the 70 /zm limts on a distant Kuiper Belt. While the uncertainty in 
this simple model is large, the total mass of comets required to maintain the hot material 
argues against the cometary hypothesis. 



4-4-4- A "Super- Comet" as a Source of Small Grains 

An alternative to a swarm of comets is a single giant comet scattered into the inner 
HD69830 system, perhaps an object like Sedna with a radius of 1,000 km and composed of 
ice and rock. Such an object would have a surface area (1000/25) 2 = 1600 times that of 
Hale-Bopp and in an low eccentricity, 0.5 AU orbit would evolve gas and small dust grains at 
a rate dM (t) / dt— 1600 times greater than Hale-Bopp. In an idealized picture, this material 
would build up over r PR , would last until the entire comet evaporated in a time, r evap , and 
then decline in t pr . During r evap the average amount of material would be equivalent to 
dM(t)/dt x t pr ~ 1.6 x 10 6 Hale-Bopps which agrees reasonably well with the amount of 
emitting material estimated in §4.4.3. 

We can estimate the evaporation time by noting that the mass loss rate is proportional to 
the surface area of the comet: dM(t)/dt = An<kR(t) 2 where $ is the flux of material emitted 
in g cm~ 2 s" 1 and R(t) is the radius as a function of time. Then, with R(t) = (?>M(t) / Aitp) 1 ^ 
and M(0) = A/3npR(0) 3 , we get 

= 4,(fl(0)p-«)» 
6p A 

and r evap = R(0)p/§. At perihelion, Hale-Bopp gave up some 200 metric tons sec" 1 in 
gas and small grains (Lisse et al. 1997) and possibly as much as 2,000 metric tons sec -1 in 
larger grains detected at submillimeter wavelengths (Jewitt & Matthews 1999). For Hale- 
Bopp with an assumed 25 km radius and a mass loss rate of 200 tons sec -1 for small grains, 
$ = 76 g cm -2 s -1 . An escape velocity from a Sedna-like object is about 0.5 km s _1 , half that 
of Pluto's, and comparable to the thermal velocity of H 2 at 300 K. Small grains entrained 
with the gas could escape the nucleus while larger grains would probably fall back under the 
self-gravity of the object. Therefore, we adopt the lower, small-grain mass loss estimate for 
this putative super-comet. It is worth noting that this grain segregation mechanism would 
account naturally for the lack of large grains with their associated long wavelength emission. 

A 1,000 km sized comet trapped in a low eccentricity orbit (making <3> roughly constant) 
would have an evaporation lifetime of r evap = 2 Myr and could provide material over a long 
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enough time to have a reasonable chance of being observed by Spitzer. Effects of self-gravity 
and ice-rock differentiation would reduce the flux of evolved material, but would extend the 
duration of the outburst. This scenario would be less attractive if the super-comet were on 
a highly eccentric orbit, spending much of its time away from the star. The requirement for 
the super-comet to travel in a low eccentricity orbit close to the star suggests the presence 
of a planet located close to the star to effect its capture. 

Finally, we note that the rupture of a single Sedna-like object in the inner HD69830 
system could release a swarm of 10 5 — 10 6 smaller, 10-25 km bodies which in turn could 
produce a large number of small grains as envisioned in the preceding section. Since the 
small grains would last only tpr, a significant fraction of the swarm would have to survive 
for a few million years to replenish the small grains and thus to make the outburst last long 
enough for us to have a chance of seeing it. 

In either case, the Sedna-like object could have been dislodged from a stable, more dis- 
tant orbit by a planet in the HD69830 system or by the effects of a passing star. The Sedna- 
like object might even have been captured from the passing star itself (Kenyon & Bromley 2004), 
although the space motion of HD69830 (Eggen 1998) indicates no nearby approaches to any 
known star within the past 20 Myr (V. Makarov, private communication, 2005). 

4-4-5. An Observational Test for the Origin of the HD69830 Excess 

There is a powerful observational test that could decide between the two models. The 
cometary scenario would bring a great deal of volatile material into the inner HD69830 
system. In Hale-Bopp, species such as CO, H 2 0, HCN, and oxygen were detected at wave- 
lengths from the millimeter to the ultraviolet. Sensitive spectroscopy of HD69830 might 
reveal evidence for volatiles seen either in emission or in absorption against the star, simi- 
lar to what has been detected toward f3 Pictoris in the visible (Lagrange-Henri et al. 1992; 
Petterson & Tobin 1999) and in the UV (Vidal-Madjar et al. 1994; Roberge et al. 2000), 
and toward HD100546 (Vieira et al. 1999; Bouwman et al. 2004). While examination of 
the IUE archive spectrum of HD69830 (Hufnagel & Smith 1994) shows no features attributable 
to gas in a shell around the star as is seen in (5 Pictoris, e.g. excess Mg II absorption at 279.5 
and 280 nm, more careful searches for volatile material are clearly warranted. Detection of 
such volatiles would strongly support the cometary hypothesis and represent a milestone in 
our understanding of how organic molecules are transported from the cold outer solar system 
into the region of the terrestrial planets. 
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5. The Detectability of Planets in the HD69830 System 

Finally, we note that the intense diffuse emission toward this star makes it unlikely that 
future astronomical telescopes will ever detect directly the light from terrestrial planets in 
the HD69830 system. For telescopes like the Terrestrial Planet Finder (Beichman 2003), the 
photon noise at either visible or mid-IR wavelengths associated with a target's exo- zodiacal 
cloud begins to dominate other noise sources at levels exceeding 10 times the brightness of 
our own zodiacal cloud. Unfortunately, the HD69830 cloud, being more than 1,000 times 
brighter than our cloud, will hide planets like our own from view. 

From the standpoint of TPF it is useful to note that the limiting IRS sensitivity described 
herein will allow us to detect or set limits to small grains located in the "habitable zones" at 
least a factor of 10 fainter than that around HD69830. A focused program of IRS observations 
should reach ~ 100 times the level of our zodiacal cloud and thus provide a good initial filter 
for TPF target selection. 

6. Conclusion 

We have identified an excess of emission between 8 and 35 /im toward the nearby K0 V 
star HD69830. The IRS spectrum of the excess reveals prominent features due to crystalline 
silicates such as forsterite, which must be present in small, sub-micron dust grains for the 
features to be so prominent. The lack of longer wavelength emission and the ratios of various 
spectral features suggest that most of this material has a temperature in excess of 400 K and 
must be located within 1 AU of the star. An excess of this type (strong at 25 /xm, but weak 
at 70 /im) is unique among the sample of 84 main sequence stars examined by Beichman et 
al. (2005) and in other IRAS/ISO samples, implying a short lived event or a low duty cycle 
process. 

Estimates of the mass in an enhanced asteroid belt needed to replenish the small dust 
grains seen by Spitzer are quite uncertain: a) 64 times the amount in our solar system based 
on an extrapolation of the emitting surface area of small grains up to 10 km bodies; or 
b) 22 times the amount in our solar system based on dust production in asteroid-asteroid 
collisions. If the mass loss rate inferred from the present day excess were constant over 2 
Gyr, the system would have lost at least 500 times the amount of asteroidal material in 
our solar system. It is more likely that the excess is sporadic and the total amount of lost 
"asteroidal" material is smaller than this value. 

Alternatively, small grains could arise from a single, Sedna-sized "super-comet" recently 
captured in the inner HD69830 system and "evaporating" over a few million years or breaking 
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apart into smaller objects. This super-comet would have to be trapped in a low eccentricity 
orbit around 1 AU from the star. Spectroscopic follow-up to search for volatiles associated 
with comets would help distinguish between these two hypotheses. 
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Fig. 1. — The Spectral Energy Distribution of HD69830 shows a Kurucz photospheric model 
(black) fitted to short wavelength data from Hipparcos and 2MASS (red). The 2MASS data 
(J, H, K a ) are quite uncertain due to the brightness of the star. Corrected spectra from the 
three IRS modules (green) are shown along with MIPS and IRAS photometry (red). 
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Fractional Excess Over Photosphere 




Wavelength (micron) 

Fig. 2. — After correction with the "super-flat" small deviations from a smooth photosphere 
are readily detectable in these IRS spectra. Shown are the fractional excesses, (F U) obs — 
F V} Kurucz) / Fv,Kumcz for 2 stars with no evidence for excess at any wavelength, HD173667 
(red) and HD142373 (green) and for HD69830 (black) which shows prominent excesses from 
8 to 35 /im. 
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Fig. 3. — Top) The spectrum of the excess of HD69830; bottom) for comparison, the spec- 
trum of the comet Hale-Bopp (Crovisier et al. 1996) normalized to 400K as described in the 
text. Identifications of some of the mineralogical features (Table 2) are indicated. 
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Fig. 4. — As described in the text, a simple model using crystalline silicate grains in a disk 
extending to ~ 1 AU (dashed blue line) is fitted to the IRS spectrum of the excess (solid 
line) . 
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Table 1: Infrared Photometry of HD69830 

Facility A(/im) F v (Jy) Photosphere(Jy) 1 Significance 2 , \ 
IRAS 12TfJ 0.68 ±0.041 fX61 L7 
MIPS 23.7 0.23 ±0.012 0.16 6.1 
IRAS 25.0 0.24 ±0.026 0.14 3.8 
MIPS 70.0 0.019±0.003 0.018 0.4 



Note. — 1 From Kurucz model fitted to short wavelength data 
Note. — 2 x ^(Observed- Kurucz) /-^/noise 2 + (model uncertainty) 2 . 
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Table 2: Features in Excess Compared with Identifications in Comet Hale-Bopp 1,2 



A(/xm) 



Comment 



Identification 



11 (9.4-12.0 ) 
10.2 

11.0-11.3 

16.7 (16.2-17.5) 
19.3 (17.7- 20.7) 

23.8 (22.7-24.7) 
27.8 

33 -35 



Very strong emission 
Shoulder within 11 um feature 
Peak within 11 um feature 
Broad feature on shoulder of 19 fxm feature 
Very broad emission feature 
Prominent, broad feature 
Weak feature or plateau 
Weak feature 



Crystalline Forsterite 
Crystalline olivine 
Crystalline olivine 
Crystalline olivine 
Crystalline olivine 
Crystalline olivine 
Crystalline pyroxene (Enstatite) 
Crystalline olivine, low SNR 



Note. - 



1,2 Crovisier et al. 1996, 1997; Wooden et al. 2000; Henning 1998. 
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Table 3: Properties of Best Fit Model 

Model 

Best Fit Model Parameters 

Density power law index, a -0.4 
(constrained) 

W (AU) 1.0 
Surface area of Hale-Bopp like grains 1 (cm 2 ) 2.7 x 10 23 

Reduced X 2 (250 d.o.f) 3^ 

Derived Parameters 

Predicted 70 /zm flux density (mJy) 0.75 
Dust Mass in 0.25 grains (M e ) 1.0 x 10~ 9 



Note. — 1 Assumes optical constants derived from scaled Halc-Bopp spectrum with amorphous silicates 
subtracted to remove the 9.3 ^m feature seen in Hale Bopp, but not in HD 69830. 
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Table 4: Comparison of Solar System and HD69830 Model 



Parameter 



Solar System 



HD69830 



System Parameters (see text) 
Orbital Distance (AU) 
Stellar Mass (M ) 
Stellar Luminosity (L Q ) 
Grain Size (/xm) 

Fractional Dust Surface Density 

(enhanced after collision, see text) 
Fractional Dust Surface Density 1 (>30 AU) 



10 



-6 



3 
1 
1 

10-30 

10~ 7 (<3 AU) 
6 x 10~ 7 
-10~ 7 (estimate) 



3 x 10" 



1 

0.8 
0.45 
0.25 
:<1 AU) 

< 6 x 10~ 6 



Timescales for Dust Destruction ( 1 AU) 
Tcoii (yr) 
Tpr (yr) 



lxlO 6 
2-6xl0 5 



400 
700 



Asteroid Belt Mass (a max = 10 km; M 0y ) 

Our Solar System 2 

From emitting surface area 3 

From mass production rate 4 

From integrating mass loss rate over 2 Gyr 



0.014 x 10" 



0.9 x 10~ 3 
0.3 x 10~ 3 
- 10 x 10~ 3 



Note. — 1 Stern (1996) for outer solar system and from 70 /im limit for HD69830; 2 From Bidstrup ct al. 
(2004) and Krasinsky et al. (2002); 3 From scaling the emitting surface area to total mass (§4.2); 4 Based on 
asteroid collisions for PR-drag dominated dust, n gr oc Af? (§4.4.2). 



